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Abstract 

We obtain exact solutions of the one-dimensional Schrodinger equation for some 
CN ■ families of associated Lame potentials with arbitrary energy through a suitable ansatz, 

which may be appropriately extended for other such a families. The formalism of 
' supersymmetric quantum mechanics is used to generate new exactly solvable potentials. 
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1 Introduction 

q-i The solution of the one-dimensional Schrodinger equation with periodic potentials is im- 
portant due to the possibility of finding interesting models which could be used in physics. 
This kind of potentials admits a spectrum composed of allowed energy bands, in which the 
physical eigenfunctions are bounded, separated by energy gaps, where the eigenf unctions are 
unbounded and thus they cannot have physical meaning due to their exponential growing 
when we move far away from a given position. 

Of special interest are the so-called exactly solvable problems which, unfortunately, in 
the periodic case include very few potentials, e.g., Lame and some others. To be precise, 
by exact solvability we mean that for the corresponding potential it is possible to determine 
analytically the physical eigenfunctions for energies in the allowed bands (the edges included), 
as well as the non-physical solutions for energies in the gaps. It is important to notice 
that in the direct spectral problems the unphysical solutions are neglected because they are 
apparently useless. However, in the inverse spectral problems, which somehow include the 
supersymmetric quantum mechanics as a particular case [1,2], they can be used as seeds to 
generate exactly solvable potentials from a given initial one [3-8]. 

In recent times it has been realized that the associated Lame potentials, which include in 
particular the Lame case, admit explicit expressions for the band edge eigenfunctions [9-13]. 
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Thus, it is natural to explore if those potentials are exactly solvable in the sense pointed 
above. If the answer turns out to be positive, then the band edge eigenfunctions as well 
as the non-physical solutions for the gaps can be used through super symmetric quantum 
mechanics to generate new exactly solvable potentials, some of which could provide some 
interesting physical information. 

In the next section we will consider some particular families of the associated Lame 
potentials and show that they belong to the exactly solvable class. Then, the formalism 
of supersymmetric quantum mechanics will be applied to generate new exactly solvable 
potentials. We will end the paper with our conclusions. 



2 Associated Lame potentials: general solutions 



The Schrodinger equation for the associated Lame potentials in Jacobi form may be expressed 
as 



d 2 ip(x) 
dx 2 



+ V(x)ijj(x) = Eijj(x) 



x G 



en x 

V{x) = m(m + l)fc 2 sn 2 x + £(£ + l)k 2 — ^~ 



dn 2 x 



(1) 



Our aim is to find exact solutions of Eq. (P) for arbitrary values of E. Here snx = 
sn(x,k), cnx = cn(x,k), dnx = dn(x,k) are three Jacobian elliptic functions of real 
modulus k 2 (0 < k 2 < 1) and of double periods AK,2iK';AK,AiK'; 2K,UK' respec- 



2 sin 2 . 



tively, K = f£ /2 d(j)/y/l — k 

K' = £ /2 d(j>/y/l - V 2 sin 2 0, k' 2 = 
step of the process we express Eq. 
from the Lame potential: 



is called complete elliptic integral of second kind and 

1 — k 2 is known as complementary modulus. In the first 
JTJ) in Weierstrass form to easily observe the difference 



d 2 jj 
~~d^ 

where z = (x - 
1). In Eq. (J2D, 



+ 



m(m + l)p(z) + 



l(l+l)e 2 e 3 
p(z) - e x 



ip = Eip, ti = ei 



2,3, 



(2) 



- iK')/yfe e C, %j){x) = ip(z(x)) and E = [E — £(£ + l)]e 3 + e 3 m(m + 
p(z) = p(z; uj, oj') is the Weierstrass elliptic function of half-periods lo = 
Kj y/Ez, uj 1 = iK' I y/e^ and p(u>i) = e«, uii = u, uj 2 = uj + uj', uj 3 = uj'. It is clear now that the 
differential equation (J2J) has regular singularities at the poles of p(z) and at points congruent 
to z — oui, the latter being absent in Lame equation which corresponds to i = 0. Our aim 
is to find two linearly independent solutions ifji(z), i = 1,2 of Eq. (J2J). It is well-known that 
*&(z) = ipi(z)ip2{z) satisfies a third order differential equation (see e.g. [14]) 



d 3 ^ 
~d£ 
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m(m + l)p(z) + 
P( 

m{m + 1) 



+ l)e 2 e 3 



p{z) - ei 
+ l)e 2 e 3 



-E 



dz 



{p(z)- ei y 



p'{z)V = 0, 



(3) 



where throughout this article the prime (except in K' and uj') will denote derivative with 
respect to the shown argument. To obtain its solutions we propose an ansatz and use a 
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fitting procedure to automatically fix m, £ and also the unknown quantities, namely; 



9(z) = [p(z)-e 1 ] + A 1 + f 2 . (4) 

p{z) - e 1 

Let us notice that this ansatz is motivated by realizing that Eq. (j3J) has an additional regular 
singularity at z = u)\ compared to Lame case {£ — 0), and so it is expected that its solution 
will also have a singularity at this point. Substituting (JH) into (J3J) and equating coefficients 
of each power of [p(z) — e x ] to zero, we find the following solutions for A x , A2, m and I: 

a) m = 1, -2; £ = 1, -2; A x = E + e u A 2 = e 2 e 3 , (5) 

b) m = 1,-2; £ = 0,-1; A x = E + ei, A 2 = 0. (6) 

The solutions (0) and © correspond to eight points in the m — i plane, four of which (for 
t = 0, —1) lead to a Lame potential. Of the other four points, the one (1,1) lies in our 
effective region (m, I > — 1/2, m > £). Thus, the ansatz (@J) will give us the general solution 
of the associated Lame potential with m = £ = 1 for arbitrary E. 

The ansatz (j3J) may be extended by inserting additional terms involving higher positive 
and negative powers of [p(z) — ei]. For instance, one may take 

¥(*) = ip(z) - ei ] 2 + B 1 [p{z) - e x ] + B 2 + ^ - , (7) 

for which the solutions become 

c)m = 2,-3; £= 1,-2; B 1 = 2e 1 + ^, B 2 = ( % - eA B x , B 3 ' 2< * Bl 




3' 3 ° 3 



E IE 

d)m = 2,-3; £ = 0,-1; 5 X = 2ei + -, B 2 = I - - e x ) B x + e 2 e 3 , 5 3 = 0. (9) 

Once again, out of the eight points four correspond to a Lame model and in our effective 
region we have the point (2, 1). 

We have found that the product of solutions of the associated Lame equation (J2J) takes 
the form: 

I 2 

■*( z ) = — l[[p(z) - p(a r )\, m = £ = l, (10) 

1 3 

*(z) = — n^)"^)]' rn = 2,£=l. (11) 

In Eqs. firU|) . (fTT|) the quantities p(a r ), r = 1,2 and p(b r ), r = 1,2,3 are the zeros of the 
polynomials arising in the numerators of and (J7J), and this implies that the determination 
of a r and b r involves transcendental equations of the type p(t) = c. Since p(z) is an 
even function, the resulting ambiguity of signs in a r , b r has to be fixed from the convention 
*lf'(z) > at the points z = a r , b r . 



3 



Once we know the product of solutions, it is straightforward to obtain two linearly inde- 
pendent solutions for the associated Lame equation ((TJ) following the same procedure adopted 
for Lame [14]. Up to some constant factors, our final results are as follows: 



1. V(x) = 2k 2 sn 2 x + 2k 



2 cn 2 x 
dn 2 x 



(12) 



n 
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-iK' 
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x-iK' 



exp 
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=F- 
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r=0 



2.V(x) =6k 2 sn 2 x + 2k [ 



, cn 2 x 



dn x 



n 



r=l 



(T 



x-iK 7 



±6 r 



x-iRT' 



+ cr 



x— iK' 



exp 



£c(V 



r=0 



(13) 

(14) 
(15) 



In above solutions we have taken a = b = ^U\, an d C( x ) are the quasi-periodic 

Weierstrass elliptic sigma and zeta functions, which are defined by C( x ) = ~p{ x )i [l n a ( x )]' = 
C(x). It is not very difficult to check that both solutions ((H)) become proportional to the 
three band edge eigenf unctions for the associated Lame potentials with m = £ = 1 when E 
takes the three eigenvalues Eq^ = 2 + k 2 =p 2 A;', — 4, which is consistent with the fact 
that this potential has one finite energy band and one finite gap [10,12]. The realization 
that the associated Lame potential for m = I — 1 and modulus parameter k 2 is exactly 
solvable was recently found by noticing [13] that it can be obtained from the Lame one 
[take m — 1, £ — and a transformed modulus parameter in Eq. (JTJ] via a well-known 
Landen transformation [14]. In contrast, here we have proved that property by directly 
finding the general solution for arbitrary Bel. On the other hand, up to constant factors 
both solutions (|T3|l reduce to the five band edge eigenfunctions for the associated Lame 
potentials with m = 2, 1=1 when E takes the five eigenvalues Eq = 4k 2 , E\^ = 5 + 
k 2 =p 2V4-3F, E 2 , 3 = 5 + 2k 2 T 2V& 4 - 5A; 2 + 4 [10, 12]. This is again consistent with the 
spectral properties for this associated Lame potential which has two finite energy bands and 
two finite gaps. To the best of our knowledge, the fact that the associated Lame potentials 
for m = 2, £ = 1 are exactly solvable was previously unknown. 



3 Supersymmetric partner potentials 

In the modern approach to the first-order supersymmetric quantum mechanics (SUSY QM), 
in which a first order differential intertwining operator is used to implement the transfor- 
mation, the seed Schrodinger solution u(x) can be either physical or non-physical but it has 
to be nodeless to avoid singularities in the new potential V(x) = V(x) — 2[lnit(x)]" (see the 
collection of articles in [2]). This is achieved by taking the factorization energy e such that 
e < E , where E is the lowest band edge energy. In particular, if we chose u(x) as any of 
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the two Bloch functions ipx^i^) derived in the previous section then V(x) will be periodic, 
with the same band spectrum as V{x). On the other hand, if we choose u(x) as a nodeless 
linear combination of the two Bloch functions associated to e, then V(x) will present 

a periodicity defect, and the spectrum of H will consist of the allowed energy bands of H 
plus an isolated bound state at e, for which the corresponding eigenfunction l/u(x) will be 
square-integrable [15-18]. 

Let us mention that for the second order SUSY QM, which involves differential intertwin- 
ing operators of second order, the key function which has to be nodeless is the Wronskian of 
the two seed solutions Ui(x), U2(x) associated to the factorization energies e±, t% [8,16-20]. 
In this frame it is possible to use solutions for which e±, 62 have unexpected positions, e.g., 
both can lie in a finite energy gap and produce however non-singular SUSY transformations. 
To be brief, in this letter we will not apply the second order SUSY transformations and we 
will restrict our discussion to the first order cases mentioned above (see however [21]). 

By applying now the first order SUSY algorithm to the associated Lame potentials with 
m — £ — 1 and m — 2, £ — 1, using as seed solution any of the Bloch functions given 
by Eqs. (|13p and (JT5j) respectively, we arrive to the following new exactly solvable periodic 
potentials which are strictly isospectral to the corresponding associated Lame potential: 



1. m = 

V(x) 



-Ak 2 sn 2 x + 2 



+2 fc 2 



r=l 
2 



snx cnx dnx ± f3(a r ^fe^) 



2 + ^a 2 (a rV ^) 



r=l 



+ 2 



(16) 



2. m 



3 r 



V(x) 



-Ak 2 sn 2 x + 2^2 



+2 U 2 



r=l 
3 



snx cnx dnx ± f3{b r ^/e^) 



sn 2 x + a 2 {b r y/e^) 
3 + J2a 2 (b r ^T 3 ) 



r=l 



+ 3 



;i7) 



where a 2 {r) = — l/(/c 2 sn 2 r), /3(r) = — cnr dnr/(A; 2 sn 3 r). 

On the other hand, if the seed solution is a linear combination of the two positive definite 
Bloch functions, u{x) = ip\(x) + A^^), then for A < 0, u{x) will have always a node, 
inducing then a singular SUSY transformation. For A = and A = 00 we will recover once 
again the previously discussed case when u{x) is one of the Bloch functions. Finally, it is 
interesting to observe that for A > 0, u(x) will be nodeless, and the new potential will not 
be strictly periodic (it will have a periodicity defect). The spectrum for the SUSY generated 
potential V(x) will consist of the allowed energy bands of the initial potential plus an isolated 
level embedded in the infinite region below the 'ground state energy' of the initial potential. 
Unfortunately, here the explicit expressions for V(x) are not compact. Thus, we decided 
to illustrate this case graphically, and an example of the SUSY partner potential with one 
periodicity defect of the associated periodic Lame potential for m = i = 1 is given in Fig. 1 
while for m = 2, £ = 1 the corresponding SUSY partner potentials are shown in Fig. 2. 
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Figure 1: Potential V{x) with a periodicity defect (black curve) generated from the associated 
Lame potential with m = I = 1 (gray curve), k 2 = 0.99, and u(x) = V>i (x) + 3^2 (s) /2 with e = 2.4, 
which leads to a\ = —1.089, a-i = 2.607. The new Hamiltonian H has an extra bound state at e. 

4 Conclusions 

Through an appropriate ansatz, we have solved the stationary Schrodinger equation for 
the associated Lame potentials with an arbitrary energy corresponding to the parameter 
pairs (1,1) and (2,1). This suggests that the associated Lame equation with any integer 
values of the potential parameters is exactly solvable. This assertion can be proved case 
by case by appropriately modifying the ansatz, given by Eqs. (JI},© and using a fitting 
procedure to find the corresponding analytic solution. For instance, one may fit the solution 
^(z) = 5^r=-2 C r [p{z) — ei] r to (JHJ), which will effectively correspond to the point (3,2) 
in the m — i plane, and so on. However, one of the aims of this paper was to show that 
this can be done through the simplest non-trivial available examples. On the other hand, 
the first order supersymmetry transformations were used to generate new exactly solvable 
potentials which can be either periodic or with a periodicity defect, depending on how we 
choose the seed Schrodinger solutions. This kind of SUSY transformations, together with 
the second order ones, represent the most simple theoretical tools for designing potentials 
with prescribed spectra, a subject which every day seems closer to physical reality. 
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Figure 2: Potential with a periodicity defect (black curve) produced from the associated Lame 
potential with m = 2, i = 1 (gray curve), k 2 = 0.95, and u(x) = ipi(x) + ^2 0*0 with e = 3.5, which 
gives b\ = —2.392, 62 = 1-26 + 0.614 i, = b\. The new Hamiltonian H has an extra bound state 
at e. 
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